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Abstract 

Graphene intercalation materials are potentially promising for the implementation of the 
ultra-low power, excitonic-condensate-based Bilayer pseudoSpin Field-Effect Transistor (BiS- 
FETs) concept, as well as other novel device concepts requiring a graphene interlayer dielec- 
tric. Using density functional theory (DFT) we study the structural and electronic properties of 
bilayer graphene intercalated with iodine monochloride (ICI) and iodine monobromide (IBr). 
We determine the structural configuration of ICI and IBr graphene intercalation compounds 
(GICs). We also conduct an in-depth exploration of inter-layer electronic coupling, using ab 
initio calculations. The presence of intercalants dopes the graphene layer. It also reduces, 
but does not eliminate, the electronic coupling between graphene layers, which may enable 
BiSFET operation. In addition, we present experimental results for IC1-GIC synthesis and 
characterization. 
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Introduction 



Graphene has generated keen interest for novel electronic applications such as the Bilayer pseu- 
doSpin Field-Effect Transistor (BiSFET). 2 The BiSFET is an ultra-low power device concept 
based on gated double-layer graphene surrounded by low-K" dielectrics to produce an excitonic 
Bose Condensate. The BiSFET requires substantial reduction of the coupling between the two 
graphene layers, but coupling should not be entirely eliminated. One potential way to achieve 
this goal is to insert intercalants between graphene bilayers, or into graphite. This insertion of 
intercalant between successive layers of graphene forms what are known as graphene intercala- 
tion compounds (GICs).- The choice of intercalant species and amount of intercalation provides 
a useful knob to control various properties like carrier concentration and interlayer electrical con- 
ductivity. Bilayer graphene (BLG) is made up of two separate layers of graphene stacked usually 
in AB and sometimes in AA configuration. Native BLG has a very strong inter-layer interaction 
and thus high c-axis conductivity making it unsuitable for the BiSFET. The presence of intercalant 
between the carbon layers makes their interaction indirect, thereby reducing its magnitude. 4 In the 
case of BiSFET, achieving optimal coupling is critical for the formation of the condensate. GICs 
are also being considered for energy storage, 5 nanoelectronics 6 and spintronics. 7 In order to de- 
sign these novel devices it is critical to understand the structural and electronic properties of these 
intercalation compounds, specially inter-layer coupling. 

In this paper we study intercalation of graphene with iodine monochloride (IC1) and iodine 
monobromide (IBr). IC1 and IBr both form acceptor stage 1 intercalants with graphene 3 and 
are very promising for our devices. The I-Cl(Br) form covalent bonds through the p electrons 8 
and graphene it orbitals are modified by interaction with iodine % orbitals. 9 There is consid- 
erable charge transfer from graphene to the halogen molecules and also from iodine to chlo- 
rine(bromine) in ICl(IBr)-GIC. The reported c-axis conductivity is very low with an anisotropy 
ratio of pel pa = 10 3 .— Sugihara 11 has proposed a theory of c-axis conduction in acceptor-type 
GICs via hopping and the low c-axis conductivity points to strong localization of % electrons in 
the graphene layer. All of this makes IC1-GIC and IBr-GIC strong candidates for the formation 
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of excitonic condensates in graphene double layers. We first address the question of structural 
configuration of these intercalation compounds. We observe that while the iodine sites proposed in 
earlier literature are correct, the chlorine and bromine sites are not. In fact there is a possibility that 
the intercalant is incommensurate with graphene. Then, electronically, we observe hole-doping of 
graphene along with a small gap opening associated with breaking the sublattice symmetry within 
the individual graphene layers, as well as a smaller energy splitting associated with the graphene 
interlayer coupling. The interlayer coupling is mediated via the the overlap between C-p and the 
intercalant orbitals. The latter represents a drastic but not complete reduction of electronic cou- 
pling between the graphene layers. These results are promising towards the realization of a GIC 
based BiSFET. 

Structural and Computational details 

Iodine monochloride (IC1) and iodine monobromide (IBr) are molecular species that form accep- 
tor stage 1 intercalants with graphene. 3,8 While being similar to Br2-GIC, 12 which has been well 
studied, the structural properties of these intercalation compounds are still subject to debate. In this 
paper we hope to illuminate that issue. In line with most of the literature on IC1-GIC, as well as our 
own experiments, we constructed models corresponding to vapor phase intercalation of graphene 
with IC1 and IBr. Heerschap et al. studied the micro structure associated with in-plane imperfec- 
tions, noting that the C planes in graphite-ICl are AA stacked. 13 The in-plane structure for graphite 
intercalated with IC1 was first reported by Tumbull et al— in 1966. They studied diffraction pat- 
terns to conclude that ICl-to-graphite layer spacing is 3.65 A while the c axis parameter is 7.0 A. 
They also suggested a monoclinic lattice for the intercalant with lattice vectors a ~ 58 A and b 
~ 28 A. In parallel they proposed a smaller pseudo-cell structure that fit their diffraction results. 
This latter structure had parameters a = 4.92A, b = 19. 2A and Y = 93.5°. While the positions of 
iodine were obtained by the Sayre-Cochran relation, 14 the chlorine sites were obtained by trial and 
error and the assumption of centra symmetry. This structure was however challenged by Moret et 
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al.— in 1983, after they performed a high resolution photographic study of IC1-GIC suggesting 
that the structure was incommensurate. A different study done by Ghosh et al— in 1984 found 
that a much bigger unit cell fit the diffraction pattern better. This unit cell was nearly double the 
size of the Turnbull model along the b lattice vector. Wortmann, Krone and Kaindl— ^ followed 
up with a model which was a refinement of the Turnbull structure. They specified that there was 
parallel alignment of intercalated IC1 molecules with respect to the graphite planes. Moreover the 
intra-molecular I-Cl bond length was 2.54 A, considerably larger than that in the gas or solid phase. 
Keeping the same iodine sites as Turnbull et al. and using their new I-Cl bond length along with 
an assumption about inversion symmetry, they proposed a new pseudo-cell which was the basis of 
our initial structure in this study. 

The density functional theory (DFT) calculations were done using the Vienna Ab-inito Sim- 
ulation Package 1 — - — with PAW pseudopotentials 21 and a Generalised Gradient Approximation 
(GGA) 22 based exchange correlation correction. Our model constitued of 36 C atoms in one C 
layer with 4 IC1 units per unit cell. Taking C-C bond length as 1.417 A we get the lattice param- 
eters as a = 19.12953c- Lilly, b = 4.908? and T = 93.670°. We took the c lattice constant as 17 
A with 10 A of vacuum between neighboring supercells. The distance between IC1 and either C 
layer was 3.5 A with IC1 bond length as 2.5544 A (Q~|). Starting with this structure we optimized 
the Kmesh to 2x8x2 and energy cut off to 450 eV while checking for convergence. Subsequently 
we relaxed the ions and optimized the in-plane lattice constants until the forces were smaller than 
0.01 eV/A. 

After relaxation, the structure obtained for IC1-GIC (JD was slightly different from the start- 
ing structure. To begin with, there was a small volume expansion such that the lattice constants 
corresponding to the relaxed structure were 1.0036 times the initial lattice constants. The relaxed 
graphene layers displayed some buckling, as is evident fromCQ The equilibrium interlayer distance 
between the two C layers is 8.2 A which is larger than the ~7 A obtained experimentally. How- 
ever, we can attribute this to the absence of Van der Waal's force in our simulations. The iodine 
sites in our relaxed structure however agree remarkably well with the Turnbull model. We also see 
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Figure 1 : Model structures for IC1-GIC and IBr-GIC (a) Top and side view of IC1-GIC model 
before relaxation (b) Top and side view of IC1-GIC structure after relaxation (c) Top and side view 
of IBr-GIC structure after relaxation. 
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evidence of inversion symmetry with 2 equivalent I and CI sites in the unit cell which is in keeping 
with the spectra studied by Tiedke et al. 9 The I-Cl bond length is around 2.5 A and more. The CI 
sites are quite different from those proposed by Turnbull et al. or by Wortmann et al. . This could 
point to the incommensurability of the IC1 structure with the graphene plane. 

The reported structure and properties for IBr intercalation compounds are described as being 
very similar to those of the IC1 intercalation compounds. 23-25 Hence we started from the same ini- 
tial structural model as used for IC1-GIC and replaced CI with Br. The Kmesh was taken as 2x6x2 
and energy cut off as 450 eV as they both sufficed for convergence of total energy. Relaxation of 
ions and optimization of lattice constants was performed using VASP with PAW pseudopotentials 
and GGA for exchange correlation approximation. The resulting structure is shown in [T] Again 
there was an expansion in the unit cell with the fractional change in lattice constants being 1 .00382. 
This is larger than that seen for IC1-GIC. As a result the I-Br bond lengths are also larger than I-Cl, 
clustering around 2.7 A. The I and Br sites are very similar to the I and CI sites, respectively, with 
the structure displaying inversion symmetry and 2 inequivalent site type for I and Br. We thus ex- 
pect the electronic and transport properties of both these intercalated compounds to be qualitatively 
similar. 



Experimental results 

We used IC1 to intercalate bulk-graphite and few layer graphene (FLG) into stage 1 ICl-GICs. X- 
ray diffraction (XRD) was used to determine the stage number of the bulk-graphite IC1-GIC. We 
extracted the c-axis spacing from the XRD peak positions, which increases from 3.35 A in graphite 
to 7.15 A in stage 1 IC1-GIC. We then intercalated FLG and bilayer graphene (BLG) flakes and 
used Raman spectroscopy to characterize the resulting BLG-GIC. We estimate a Fermi level shift 
of ~ 0.7 eV for the BLG-GIC. This agrees very well with our DFT results. 

Bulk-graphite IC1-GIC was prepared from natural graphite flakes (NGS Naturgraphit, 10 — 20 
mm flakes) and IC1 (Sigma-Aldrich, reagent grade, 98%) using a two-zone vapor-phase intercala- 
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Figure 2: Experimental results, (a) 0—20 XRD pattern of bulk-graphite, showing characteristic 
peaks that correspond to a c-axis spacing (dc-c) of 3.35 A. The unlabeled peaks come from the 
Al substrate holder, (b) The XRD pattern of bulk-graphite stage 1 IC1-GIC shows additional peaks 
from the (001), (003) and (005) planes. The interlayer spacing is extracted as dc-ici-c = 7.15 A. 
(c) Optical micrographs of an FLG flake before (left) and after (right) intercalation show changes 
in optical contrast. The scale bars are 50 /im. (d) Raman spectra of BLG and BLG-GIC show 
an upshift of the G-peak and a change in the 2D-peak shape towards a single Lorentzian after 
intercalation. The absence of D-peaks indicates defect free graphene. 



7 



tion process.- A chunk of bulk-graphite (~ 0.1 g) and IC1 (~ 0.5 g) were introduced into separate 
chambers of a two-chamber Schlenk tube which was then pumped for 10 minutes at a base pressure 
of ~ 3 x 10~ 3 Torr to remove any traces of water vapor from the IC1. The tube was then sealed 
and inserted in an oven maintained at 35°C for a period of 48 hours, to obtain stage 1 IC1-GIC— 
The glass tube was cooled in a refrigerator at 4°C for 30 minutes to condense the IC1 vapors prior 
to removing the bulk-graphite IC1-GIC. A higher intercalation temperature and/or a shorter inter- 
calation time were found to result in GICs of higher stage number and incomplete staging. XRD 
analysis was performed on the resulting IC1-GIC using high-speed Bragg-Brentano optics (CuK a ) 
on a PANalytical X'Pert Pro MRD system. A symmetric 6—26 scan was performed to determine 
the stage number and the interplanar c-axis spacing using Bragg's law ??. 

2d m sin6=nX (1) 

[2ta),(b) show the 6—26 XRD scans for bulk-graphite and bulk-graphite IC1-GIC, respectively. 
A c-axis interlayer spacing (dc-c) of 3.35 A is extracted for bulk-graphite by fitting the calculated 
peak positions from the (002) and (004) planes to their observed values. The dc-c value agrees 
well with prior reports from literature.—" 2 ^ The XRD scan of bulk-graphite IC1-GIC shows addi- 
tional peaks between the main bulk-graphite peaks, corresponding to reflections from the (001), 
(003) and (005) planes, which is a clear signature of stage 1 GIC. 3,27 The c-axis interlayer spacing 
(C/IC1/C) is extracted to be dc-ici-c = 7.15 A, which has been reported earlier in literature.— 
The broadening of the intercalant peaks in |2tb) suggests that even though the GIC predominantly 
consists of a stage 1 superlattice structure, it also includes minority crystal phases of graphite and 
possible higher stage GIC subdomains, possibly due to incomplete intercalation. 28 We believe 
that an optimization of the intercalation recipe would result in a more homogeneous stage 1 GIC 
formation. 

To further characterize stage 1 ICl-GICs, we intercalated FLG and BLG flakes exfoliated onto 
285 nm SiCVSi substrates. The same recipe used for bulk- graphite was used to intercalate the FLG 
and BLG flakes. (2tc) shows optical micrographs of an FLG flake, before and after IC1 intercala- 
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tion. Thicker FLG flakes show a marked change in optical contrast, whereas contrast changes in 
thinner FLG and BLG flakes are negligible. BLG flakes, however, do get intercalated, as indicated 
by changes in Raman signatures after intercalation. |2td) shows the Raman spectra of a BLG flake 
measured at 532 nm, before and after intercalation. There is a large upshift of the G-peak from 
1586cm -1 to 1623cm~ 1 (+37cm~ 1 ), which is indicative of intercalation. 29 Such a large upshift of 
the G-peak cannot be due to charge-transfer doping from IC1 adsorption alone, but must be due 
to intercalation of IC1 into the BLG flake.— Using work done by Das et al.,— we estimate the 
Fermi level shift of the BLG-GIC to be ~ 0.7 eV. This agrees very well with our theoretical pre- 
dictions using DFT The 2D-peak shape of BLG also changes significantly after IC1 intercalation. 
Pristine BLG is AB stacked and exhibits a broad 2D-peak composed of four Lorentzian compo- 
nent peaks. 32 However, after intercalation, the 2D-peak shape changes significantly, towards a 
single Lorentzian. This signifies electronic decoupling of the graphene layers into an AA stack- 
ing configuration. 33 The faint shoulder on the 2D-peak at ~ 2740cm -1 could possibly be from 
un-intercalated islands in the BLG-GIC flake. 33 In addition, there is no D-peak observed for the 
BLG-GIC, indicating defect-free graphene after intercalation. Raman spectra of FLG-GIC show a 
similar behavior, with a large upshift of the G-peak and a change in the 2D-peak shape towards a 
single Lorentzian. 

Electronic structure 

We calculated the bandstructure of relaxed bilayer graphene intercalated with IC1 and IBr. As seen 
in [3l graphene maintains its linear spectrum except very near the K point. However, it becomes 
hole doped due to the presence of the acceptor type intercalant. Moreover, for both structures 
there is a small gap between the (nominally) conduction and valence bands at the Dirac point, 
i.e. between the upper and lower cones, equal to 12 meV for IC1-GIC and 11 meV for IBr-GIC, 
which we label E g \ . There is still smaller splitting throughout the conduction band and, separately, 
throughout the valence band, the value of which at the K point we label E g 2. The value of E g \ and, 
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for each band, that of E g 2 are provided inQ] 




Figure 3: Electronic structure and gap opening at the K point for, (a) and (b), relaxed IC1-GIC, 
and, (c) and (d), relaxed IBr-GIC. The zero energy reference is the Fermi level. The M, K and K' 
points correspond to those of the Brillouin zone of the graphene layers. 

There are two sources of degeneracy breaking here, that of coupling between the graphene lay- 
ers and that of symmetry breaking between the sublattices of the individual layers due to coupling 
to the intercalant. To help isolate the effects of each, we performed band structure calculations in 
structures with one of the two graphene layers removed, i.e. for just C-IC1 and and C-IBr as shown 
in Ufa) and BJc), respectively. The results, with only the valence-to-conduction band splitting, 
show coupling to the intercalant to be the source of E g \ by its presence, and weakened interlay er 
coupling between the graphene layers through the intercalant to be the source of the E g 2 splitting 
by its absence. The coupling between the graphene layers is of particular importance for the pro- 
posed BiSFET. It is well known that the two C layers in AB stacked bilayer graphene are highly 
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coupled, which leads to their parabolic band structures near the K points and diminishes the sym- 
metry between valence and conduction bands 34 close to the Fermi level. In AA stacked bilayer 
graphene this C-C coupling leads to the two Dirac points split above and below the Fermi level. 
The presence of intercalant molecules makes the C-C interaction indirect, i.e. mediated through 
overlap with the intercalant molecules. 4 This significantly reduces the magnitude of this interca- 
tion as compared to the bilayer case. This reduced coupling, mediated through hopping between 
C-p and intercalant orbitals, in principle determines the critical (maximum) interlayer condensate 
transport in the BiSFET. We also note that band structure obtained for methane-intercalated bi- 
layer graphene 35 for both AA and AB stacking turn out to be very similar. Our results for AA 
stacking can thus be qualitatively extended to intercalated AB stacked graphene, probably due to 
indirect interaction between the C layers. However, as previously noted, the interlayer distances 
we obtained during DFT-based relaxation of the intercalated structure were obtained without Van 
der Waal's force. Because the electronic coupling between the layers is strongly dependent on 
the interlayer distance, we also conducted a set of calculations with the intercalated carbon layers 
"squeezed" down from the "relaxed" 8.2 A to the experimental value of 7.15 A and with half the 
latter between the C and the intercalant layers when considering just one graphene layer. These 
latter results are summarized in[TJ and the band structures for the intercalated structures are shown 
in|Sb) and (c). We found, as qualitatively expected, that forcing the layers closer both increased 
the interactions leading to a higher E g \ and E g 2 and somewhat increased the doping of the C layers. 

Table 1: Gap openings between intercalated graphene energy bands due to interaction of the 
graphene layers with the intercalant and with each other. E g \ is the separation between the (nomi- 
nally) conduction and valence bands; the two values of E g 2 are the splitting of the conduction band 
and of the valence band, respectively, at the K point. All gaps are in meV. "Relaxed" structures are 
obtained upon DFT relaxation absent Van der Waal's force; "squeezed" structures have the layer 
separations based on the experimental intercalated graphene layer distance artificially of 7. 15 A. 
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Figure 4: Gap opening at the K point for (a) relaxed C-IC1 (b) IC1-GIC with the experimental 
C-C distance (d C -c) of 7.15A (c) relaxed C-IBr (d) IBr-GIC with the experimental C-C distance 
of 7. 15 A. The zero energy reference is the Fermi level. The K point corresponds to that of the 
Brillouin zone of the graphene layer(s). 
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Conclusion 



In this paper we have presented first-principles density functional theory based calculations of bi- 
layer graphene intercalated with iodine monochloride and iodine monobromide forming stage 1 
acceptor-type intercalation materials. These systems are highly promising for the proposed untra- 
low power BiSFET, 2 as well as perhaps other novel devices. 5 " 7 We have studied the structural 
configurations of these materials which have been subject to debate. Moreover, we have explored 
the electronic properties of IC1-GIC and IBr-GIC, specially focussing on the inter-layer interactions 
which are highly critical for novel device functionality. In addition, we have presented experimen- 
tal results on intercalation of bilayer graphene (BLG) and few layer graphene (FLG) with IC1. 
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